1 



o 

CD 

Q 



Can baryogenesis occur on the lattice? 
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We examine the question of how baryogenesis can occur in lattice models of the Standard Model where there 
is a global 17(1) symmetry which is accompanied by an exactly conserved fermion number. We demonstrate that 
fermion creation and annihilation can occur in these models despite this exact fermion number conservation, by 
explicitly computing the spectral flow of the hamiltonian in the two dimensional U(l) axial model with Wilson 
fermions. For comparison we also study the closely related Schwinger model where a similar mechanism gives rise 
to anomalous particle creation and annihilation. 
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1. INTRODUCTION 

It is well known that fermion number is not 
conserved in the Standard Model (StM) due to 
the anomaly in the baryon and lepton currents 
which arises from the rich vacuum structure of 
non-abelian gauge theories. Topologically dis- 
tinct vacua which are characterized by their in- 
teger Chern class n are separated by a poten- 
tial barrier whose minimal height is given by 
the so-called sphaleron energy, iJgph- Transi- 
tions from one vacuum to another are accom- 
panied by a change in baryon and lepton num- 
bers, AQb,l = Nc^n, with Ng the number 
of generations. It was shown in ref. [Q that 
the tunneling transition rate at zero temperature 
is enormously suppressed by the small exponen- 
tial factor e"**^/""' « e~^^°; aw being the weak 
fine structure constant. As was first pointed out 
in ref. |^ though, the transition rate gets sub- 
stantially enhanced at large temperatures, since 
the system can travel classically from one vac- 
uum to another by thermal excitation instead 
of having to penetrate the barrier by tunnel- 
ing. This happens when the thermal energy ap- 
proaches the sphaleron energy and the transition 
rate is then proportional to the Boltzmann factor 
exp(-£;,ph(r)/T). 

The existing semiclassical calculations of the 
transition rate are based on certain assump- 



tions and it is important to check the results 
by a calculation from the first principles and 
the only possibility at present to address such a 
non-perturbative problem is to use lattice Monte 
Carlo calculations. Some first attempts in this di- 
rection have been performed in the SU(2) gauge- 
Higgs sector of the StM and in the two dimen- 
sional abelian Higgs model [Q. However these 
simulations must be extended to the full StM in- 
cluding fermions. A potential problem however 
is that a satisfactory formulation of the StM is 
still lacking. On the lattice each fermion is ac- 
companied by spurious species doublers of oppo- 
site chirality which render chiral gauge theories 
vector-like. In the past three years most mod- 
els, which have been proposed so far to overcome 
this phenomenon, have been shown to fail. The 
situation looks even worse with regard to the is- 
sue of fermion non-conservation. It has been ar- 
gued in ref. (s) that currently existing models 
are unable to reproduce the fermion number non- 
conservation of the StM, since the action which 
in general is of the form 
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and also the measure in the path integral are 
invariant under a global U(l) symmetry tpx 
exp{ia)'ipx, ipx ~* ^2;6xp(— ia). One therefore 
expects this exact global U(l) symmetry to be ac- 
companied by an exactly conserved fermion num- 
ber, and hence no baryogenesis on the lattice 
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when fermions are included in the above simu- 
lations. 

We shaU outhne in this contribution how 
anomalous fermion creation and annihilation can 
occur on the lattice despite the exact fermion 
number conservation. For the sake of simplic- 
ity we have illustrated our reasoning for the case 
of a two dimensional axial model with Wilson 
fermions, however we see no reason that the re- 
sults would not generalize to four dimensions, or 
to models with Higgs fields. Most of this work has 
been previously presented in ref. Q using stag- 
gered fermions. 

2. AXIAL AND VECTOR QED2 

The massless two-dimensional axial model is 
given by the following (real time) continuum ac- 
tion 
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S = - I Sx^F^.F^" 



Sf , 



(2) 



(3) 



where 7"'" = 71 = (Ti, 7" = —70 = —i<^2, 
75 = 0-3. We take to be an external gauge 
field; only the fermion fields are quantized. The 
model is invariant under local axial gauge trans- 
formations: ^{x) — + exp{iuj{x)^5)Tp{x) , ip^x) —^ 
i>{x)exp{iuj{x)j5), A^{x) ^_Af_,{x) + df_,oj{x), 
with gauge current — iipj^^y^Tp. The ac- 
tion is furthermore invariant under the global 
U(l) symmetry 'ip{x) — > e:xjp(ia)ip{x), tpi^) ~^ 
-0(2;) exp(— ia), with a corresponding globally 
conserved vector current j'^ = iip^^ip. Accord- 
ing to standard lore the (local) gauge current 
has to be anomaly free whereas the vector cur- 
rent becomes anomalous, df^j^ = —2q where 
q = -^e'^" Ffj^i, = dfj^C^ denotes the topological 
charge density and C^(a;) = (l/27r)e^''A,,(a;) is 
the Chern-Simons current. The axial and vec- 
tor charges are defined as Q5 — J dx^j^{x) and 
Q — J (ix^j°(x), and the Chern-Simons number 
as C = / dx^C'^{x). From the current divergence 
equation one sees that 

Q{t) - g(0) - -2{C{t) - C(0)) (4) 

which relates a time dependent change in the 
Chern-Simons number to a time dependent 



change in the fermion number. Thus a change in 
the topological charge of the gauge field AC = — 1 
gives rise to the creation of a left and right-handed 
fermion. The aim of this contribution is to inves- 
tigate whether, in lattice models which have exact 
global U(l) symmetry (i.e. the vector current is 
divergence free in contrast to the above contin- 
uum example), a change in C due to a sphaleron 
transition can still give rise to a change in the 
fermion number in accordance with eq. (^. 

A nice feature of the two-dimensional axial 
model is that it is equivalent to the massless 
Schwinger model (QED2) which is well studied. 
By performing a charge conjugation transforma- 
tion on only the right (or left)-handed fermion 
fields in eq. (||), i.e. Vfl = (V'ijC)^, ipn = 

— (C'l'V'fl)"^, 'ipL = V'l: V'l = V'l O"^^ finds that 
eq. (§) turns into 



S'f = - I d'xij'r[d^ 



(5) 



which is the action of vector QED with mass- 
less fermions in two dimensions. We put a prime 
in the sequel to label quantities of the vector 
model. Under the same charge conjugation it can 
be shown that the axial current turns into the vec- 
tor current, jg — > —iip = —j'^, which is now 
the divergence free gauge current, d^j^ = 0. The 
global vector current becomes the global axial 
current — s- — ii/i 7^75V'' = ~3i^ and is anoma- 
lous, df^j'^ — 2q(x). Correspondingly Q —Q5, 
Q5 ^ ^Q' for the charges. From the divergence 
relation for jg we can derive the following formula 

Q'S)~Q'M = 'i{C{t)-cm ■ (6) 

A change in Cern-Simons number AC = — 1 gives 
rise to the creation of a left-handed fermion and a 
right-handed antifermion. It is therefore instruc- 
tive to study the spectral flow of the vector model 
in conjunction to the axial model. 

In our earlier publication we used staggered 
fermions to transcribe the two models to the lat- 
tice 1^ . Here we shall use instead Wilson fermions 
which may be somewhat more transparent. The 
euclidean lattice path integral for external gauge 
field configuration is then given by 



j {\{V^^V^:^) eM-{SF + Sw)], 



(7) 
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with 

fix 

for the vector QED2 and 

fJ,X 

-^^,+^-f,{u;,PL + u^^Pr)^.) . (9) 

for the axial model. We use lattice units in which 
the lattice distance a — 1. U^x = exp(— i^^j.) 
is the lattice gauge field and Pr^l = 5(1 i Ts) 
are the chiral projectors. We have added in (|^) a 
Wilson term of the form 

fJ,X 

+^'.u;j''.+i.}) ■ (10) 

(and for the axial model with ip'^ —>■ ipx) which for 
external gauge fields leaves the physical fermion 
at momentum (0, 0) massless and removes the 
three species doublers, that are located in the two 
dimensional Brillouin zone at momenta (0, tt), 
(tt, 0) and (tt, tt), from the spectrum by giving 
them a mass of the order of the cutoff. 

The vector model is invariant under local vector 
gauge transformations tp^ — > ip^^^^, tpx ^x'tpx, 
U^x — > ^xU^x^x-i-fi- The Wilson term however 
breaks the global chiral invariance but it has been 
shown ref. [Q that the violation of this global 
symmetry turns, in the scaling region, into the 
usual chiral anomaly. We investigate in the fol- 
lowing how this manifests itself in the spectral 
flow of the hamiltonian when the gauge field is 
slowly changed from one vacuum configuration to 
a topologically distinct one. An earlier analysis 
of the spectral flow in vector QED2 with Wilson 
fermions has been presented in ref. 

The path integral of the axial model is invari- 
ant under the unwanted global U(l) symmetry. 
There however the (local) chiral gauge invariance 
is broken by the Wilson term. It has been shown 
in ref. ||^ that gauge invariance can be restored in 
the perturbative regime by adding a mass coun- 
terterm for the gauge field to the action, however 
a complete chirally gauge invariant lattice regula- 
tor is still lacking and is a fundamental problem of 
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Figure 1. Spectral flow for vector QED2. 



lattice gauge theory. As mentioned earlier most 
of the very promising proposals for a lattice for- 
mulation of chiral gauge theories have been shown 
to fail. Much work remains to be done in this di- 
rection, and our present study of the subtleties 
of anomalous fermion production seems to be a 
necessary step, which is also of interest in its own 
right. 

3. INVESTIGATION OF THE SPEC- 
TRAL FLOW 

In this section we study the spectrum of the 
fermionic hamiltonian in an external gauge field, 
which starts out as a vacuum gauge field, goes 
through a sphaleron like configuration and ends 
up as a different vacuum gauge field configura- 
tion that is related to the first one by a topolog- 
ically non-trivial ('large') gauge transformation. 
We shall only consider gauge fields which vary 
very slowly in time so that we can make use of 
the adiabatic theorem and determine the spec- 
trum of the hamiltonian by continuity. From the 
euclidean path integral one can derive the transfer 
operator T, Z = TrT^, where N is the number 
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of time slices in the lattice. The hamiltonian is 
defined in terms of T by T = exp{—H). We shall 
consider in the following the simple gauge field 
Ai^ = A{t), Ao^ = (i.e. Uq = 1). Then the 
time dependence of the hamiltonian is expressed 
as A dependence. The Chern-Simons number is 
given by C = —LA/2tt. With this choice of the 
gauge field configuration we can diagonalize the 
transfer operator explicitly and then determine 
from its eigenvalues the spectrum of the associ- 
ated hamiltonian. We shall furthermore use anti- 
periodic boundary conditions for the fermions in 
space which implies that the spatial momenta are 
given by p = {n - \)2ii/L, n = + 1, . . . , | 
with L denoting the extent of the lattice in space 
direction. 

3.1. Spectral flow in vector QED2 

In the vector model we derived for the enery 
spectrum e(p, A) of the hamiltonian the following 
expression 



t{p,A) = eo{p-A) 
1 



cosheo(p) = 



4sin2f 



1 



2sin2f 



(11) 
(12) 



where eo(p) are energies for free Wilson fermions. 
For momenta p and gauge fields A which are small 
compared to the cutoff tt (|i^) reduces to the con- 
tinuum expression e(p, A) ~ ±{p~ A) for linearly 
spaced modes p. The physical vacuum (Dirac sea) 
is the state of lowest energy in which all levels 
with negative energy are occupied. The spectrum 
of the hamiltonian has been plotted in fig. 1 as a 
function of LA/2tt = —C using a lattice of size 
L = 16. For clarity we have displayed only the 
states in and close to the surface of the Dirac sea. 
The L states move upwards whereas the R state 
move downwards. Midway between A — and 
A = 2tt/L the Dirac sea loses an L state and 
gains an R state (coming from the Dirac sky). 

Consider now the flow of the physical state 
1^*, A) which at A = starts out as the vacuum 
state |0, A). The occupied levels of that state are 
represented by the solid and dashed lines. Now 
assign the quantum numbers Q' and Qg to this 
state. The axial charge Q'^ is given by the num- 
ber of R states minus the number of L states 
and the fermion number Q' is given by the sum 



of L and R states minus half the total number 
of states. From this it is clear that the charges 
(5/5,* and Q'^ of the physical state are both zero 
in the initial state at A = and the final state 
at A = 27r/L. The crucial observation however is 
that fermions are excitations relative to the vac- 
uum, which implies that the change in the axial 
charge is given by 

AQ^ = (g'5 ,* '35,o)final (Qs,* '95jo)initial • 

(13) 

The axial charge of the ground state Qg q changes 
by two units when changing C from to —1 since 
it looses an L state and gains an R state, i.e. 

(Q5,o)initial = 0, but (Q5^o)final = -2. The rCSult 

is that AQg = —2 which is in concordance with 
eq. (^ since AC = — 1 in our case. Similarly it 
can be shown also that AQ' — 0. Vector QED2 



e(p,A) 
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Figure 2. Spectral flow for axial QED2. 



is gauge invariant and the vacua at A = and 
A = 2-k/L are related by a large gauge trans- 
formation. In Hilbert state, these gauge trans- 
formations induce unitary transformations on the 
hamiltonian which implies that its spectra coin- 
cide at A — Q and A = 2tt/L. We also mention 
that there is a peculiarity about the lowest state 
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in the spectrum which corresponds to the largest 
momentum. It starts out as an L state (moving 
downwards) and ends up as an R state (moving 
upwards). This change of chirahty appears to be 
possible because the Wilson term mixes L and 
R states. The chirality flip occurs however only 
for the largest momentum and hence has to be 
regarded as a lattice artefact. 

3.2. Spectral flow in axial QED2 

The energy spectrum in the axial model is given 
by the formula 



e(p,v4) 
cosh""'^ 



In- 



1 + sin A cos p + 2 sin^ | 
1 — sin A cosp + 2 sin^ | 
1 + 4sin2 f - isin^ A 



^{1 + 2 sin" 1)2 - sin2 A{1 - 2 sin^ f ) 



(14) 



where we have dropped the gauge fields from the 
Wilson term in eq. (10). Both the Wilson terms 
(with and without U fields) break the (local) axial 
gauge invariance. Eq. (|lj) reduces for momenta p 
and gauge fields A, which are small compared to 
the cutoff, to the continuum expression e(p, A) = 
±p + A for linearly spaced modes p. The energies 
e{p, A) have been displayed in fig. 2, again plotted 
as a function of LA/2tt. Each state is doubly 
degenerate since it involves an R and an L state. 
Consider now the flow of the state jvE',^) which 
starts out at A = as the vacuum state |0,A). 
Like in the vector model, the state l^*,^) does 
not change quantum numbers, i.e. (<3>i')imtiai ~ 
0, (Q*)finai — 0. The vacuum however changes 
fermion number since it looses two states half way 
between A = and A = 2Tr/L where an L and an 
R state cross the surface of the Dirac sea. For the 
vacuum state we find ((5o)initiai = 0, (Qo)finai = 
—2 and hence 

AQ = (Q^ — (5o)final ~ (Q* ~ Qo)initial ■ 

= +2 . (15) 

which agrees with the continuum result (^). In 
contrast to the vector case the spectra at A — 
and A = 2tt/L are not identical which is due to 
lack of gauge invariance. For modes near the 
maximal negative or positive momenta we see 



clearly regularization effects, while the energy of 
modes near the surface of the Dirac sea return 
roughly to their original energies. Fig. 2 shows 
that the slope of the lines decreases when increas- 
ing \p\ and finally even becomes negative. We 
should stress that such an effect has not been 
observed with staggered fermions showing that 
the regularization effects are smaller for staggered 
fermions. 

Summary: We have shown based on the inves- 
tigation of the spectral flow in the external field 
approximation, that anomalous fermion creation 
and annihilation can occur, even when fermion 
number is conserved. The crucial observation is 
that fermions are excitations relative to the vac- 
uum. The global U(l) symmetry prohibits a state 
from changing its fermion number, however noth- 
ing prevents the ground state from doing so. 
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